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ABSTRACT. 

We show how the sum of a hard singularity, F2h{x,Qq) ~ x an d 
a soft Pomeron F2p(x,Qq) ~ Const, for the singlet piece of the structure 
function F2S = F2H + Fq.p for Qq ~ a few GeV 2 , plus a saturating expression 
for the strong coupling, a s (Q 2 ) = 4tt//?o log[(Q 2 + A 2 )/ A 2 ] give an excellent 
description of experiment i) For small Q 2 , ^ Q 2 < 8.5 GeV 2 , and ii) For 
large Q 2 , 10 < Q 2 < 1500 GeV 2 if evolved with QCD. The x range is 6 x 
10 -6 ^ x < 0.04. The description for low Q 2 implies self-consistent values for 
the parameters in the exponents of x both for singlet and nonsinglet. One has 
to have a p (0) = 0.48 and A = 0.470 [ap(0) = 1.470], in uncanny agreement 
with other determinations of these parameters, and in particular the results 
of the large Q 2 fits. The fit to data is so good that we may look for signals of 
a "triple Pomeron" vertex, for which some evidence is found. 
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§1. Introduction. The HERA collaborations' 1 ' 2,3 ! have produced in the last years data with increasing 
precision for the ep deep inelastic scattering (DIS) structure function F 2 (x 1 Q 2 ) at small x, and in a wide 
range of momenta, from Q 2 ~ (including real photoproduction at Q 2 = 0) to thousands of GeV 2 . This 
allows us to study in detail not only the low x behaviour at momenta sufficiently large 1 that perturbative 
QCD be applicable, which was done in ref. 5 in a full NLO (next to leading order) calculation, but also 
enables us to address the important issue of the connection between the perturbative regime (Q 2 > 10 GeV 2 , 
say) and the region (Q 2 < 10 GeV 2 ) where nonperturbative effects are determinant. 

In ref. 5, where we send for more details, we studied fits of the large momentum region with the 
hard singularity hypothesis, i.e., under the assumption that the singlet component of F 2 , F$, behaves as a 
negative power of x as x — > so that, as discussed first in ref. 6, QCD implies the detailed behaviour 

F s {x,Q 2 ) ~ B s [a s (Q 2 )} d ^ 1+ ^x- x (1) 

x— >0 

to leading order (LO). The soft Pomeron dominance hypothesis was also considered. Here one assumes ^ 
that at a given Q 2 , <~ 1 GeV 2 one has a standard Pomeron (constant) behaviour, and then evolves with QCD 
so that, to LO, 



TP ~ C ° 



9£t 



with 



4 7 r 2 (33-2n / ) 
£= |logx|, r = log 



cxp j\/cZo£ t - dirj , (2) 



a s (Q 2 Y 



In Eqs.(l), (2) the d + , d\, do are known quantities related to the singlet anomalous dimension. The full 
NLO corrections to (1), (2) may be found in ref. 5. There it was shown that both (1) and (2) could be 
used to give reasonable fits to the large Q 2 data, with a x 2 /d.o.f. of the order of unity if fitting data with 
x < 10~ 2 or x 2 /d.o.f. <~ 2 if including points up to x = 3.2 x 10~ 2 . The hard singularity expression gave 
more satisfactory results for small x, but tended to disagree with experiment for "large" values x > 10~ 2 . 

Because of this feature it was suggested in ref. 5 that one could perhaps conclude that F$ should 
contain both a hard and a soft piece: since both (1) and (2) solve the QCD evolution equations, so will the 
sum. In the strict x —> limit, (1) dominates over (2); but if the coefficients are such that Co S> Bs then 
there will exist regions of relatively large values of x where both contributions will be comparable. In fact, 
it was pointed out that a situation similar to that occurs for real photoproduction (Q 2 = 0) where it has 
been shown! 9 ' 10 1 that an expression 2 for the cross section 

ct 7( q2 =0)p (s) = C H s x + C P + C pS - p , (3) 

provides an excellent fit to the data, up to the very high HERA energies' 4 !; and formulas similar to (3) are 
also known to work well in the small Q 2 region! 11 '. 

In the present note we pursue this line in two directions. First, we write QCD inspired phenomeno- 
logical formulas that permit us to interpolate between the large Q 2 regime, with F 2 given by a combination 
of (1) and (2) plus a nonsinglet contribution, and the low Q 2 region. Here we will show that a natural con- 
sistency condition as Q 2 — > allows us to determine the constants A = Ao, p — pa finding values remarkably 
close to those given by the fits at large Q 2 : Ao — 0.470, po ~ 0.522. What is more, the formulas provide 
excellent fits to data, with X 2 /d.o.f. =j§^ for Q 2 < 8.5 GeV 2 . 

Secondly, using these values Ao, pa, and profiting from the fact proved in refs. 6, 8 that these 
quantities are independent of Q 2 , we fix them to the values found at low Q 2 and produce a fit to high 



1 We will consider the region Q 2 > 10 GeV 2 to be that where perturbative QCD is to be trusted. This is because, at 
x — » 0, the NLO corrections are very large. 

2 We have added in Eq. (3) the contribution of the p — A2 Regge trajectory, which corresponds in DIS language to 
the nonsinglet piece of the structure function; see below. 
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momenta data with an expression 3 



F 2 = (e*)[F s + F NS ]\ 
F s (x,Q 2 ) = B s [a s (Q 2 )}- d ^ 1+ ^x- x 

4 7 r 2 (33-2n f ) 
F NS (x,Q 2 ) = B N s[a s (Q 2 )]~ dNs{1 ~ Po) x Po . 



( i 1 (4) 

exp | V^oC t - diT> ; 



[rfjvs(l — Po) is the nonsinglet anomalous dimension, a known function of 1 — po]. This fit improves substan- 
tially those obtained with (1) or (2) separately, with the addition of a single new parameter, since we took 
An, po given by the low energy calculation. We arc therefore able to give a consistent description of data, 
from Q 2 = to thousands of GeV 2 , and from x ~ 10~ 6 to x = 0.032 with, as free parameters, only the 
constants c , B Sl B N $ and Q 2 ,, plus a low energy one (see below). 

§2. The low Q 2 region. The expression for the virtual photon scattering cross section in terms of the 
structure function F 2 is 

47TCK 

CT 7(Q 2 =0)p (s) = -Q2-F 2 (x, Q 2 ), with s = Q 2 /x. (5) 

In the low energy region we should, as discussed, take the soft-Pomeron dominated expression (2) to be given 
by an ordinary Pomeron, i.e., behaving as a constant for x — > (or equivalcntly, s — > oo). So the expression 
that will, when evolved to large Q 2 yield (4) is 

F 2 = (e 2 ){Bs[a s (Q 2 )}- d +( 1+x ^x- x ° 

, (6) 
+C + B NS [a s (Q 2 )]- dNS ^^x po j. 

Note that C ^ en as the gluon component also intervenes in the evolution; the relation between C, en may 
be found in refs. 8, 5, and later in the present work. 

On comparing (5) and (6) we see that, as noted in ref. 10, we have problems if we want to extend 
(6) to very small Q 2 . First of all, 

47T 

(3 lojQ 2 /A 2 

diverges when Q 2 ~ A 2 . Secondly, Eq. (5) contains the factor Q 2 in the denominator so the cross section 
blows up as Q 2 — > unless F 2 were to develop a zero there. 

It turns out that there is a simple way to solve both difficulties at the same time. It has been 
conjectured! 12 ' that the expression (7) for a s should be modified for values of Q 2 near A 2 in such a way that 
it saturates, producing in particular a finite value for Q 2 <~ yl 2 . To be precise, one alters (7) according to 

47T 

asiQ2) /3 log(Q 2 + Af 2 )M 2 ' 

where M is a typical hadronic mass, M ~ m p ~ A(rif — 2) . . .; the value M = 0.96 GeV has been suggested 
on the basis of lattice calculations. Here we will simply set M = A = A e g, to avoid a proliferation of 
parameters. For the Pomeron term [the constant in Eq. (6)] we merely replace C — > Q 2 /(Q 2 + A 2 S ). The 
expression we will use for low Q 2 is thus, 

F 2 = {e 2 q ){Bs[a s (Q 2 T d ^+»Q- 2X s x 

O 2 ( 8a ) 
+ C Q 2 + A 2 + B N s[a s (Q 2 )r dNs{1 - p) Q 2p s- p }, 

where 

& - (Q) - ftWT^PS <8b) 



We of course include NLO corrections to (4) in the actual fits. 
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and we have changed variables, (Q 2 , x) — > (Q 2 , s = Q 2 /x). 

We have still not solved our problems: given Eq. (5) it is clear that a finite cross section for Q 2 — > 
will only be obtained if the powers of Q 2 match exactly. This is automatic by construction for the Pomeron 
term, but for the hard singlet and the nonsinglet piece it will only occur if we have consistency conditions 
satisfied. With the expression given in (8b) for a s it diverges as Const. /Q 2 when Q 2 —> 0: so we only get a 
matching of zeros and divergences if A = A , p = po with 



d+(l + A ) = 1 + A , d NS (l - po) = 1 - po- 



(9) 



The solution to these expressions depends very little on the number of flavours; for nj — 2, probably the 
best choice at the values of Q 2 we will be working with, one finds Ao = 0.470, po = 0.522. The second is in 
uncanny agreement with the value obtained with either a Regge analysis in hadron scattering processes, or by 
fitting structure functions in DIS. The first is slightly larger than the value obtained in fits to DIS^ 5 ' 6 ' which 
give A = 0.32 to 0.38. Nevertheless, these figures were obtained in fits omitting the soft Pomeron-dominated 
piece; when adding it the favoured value (see below) is A = 0.43, and, as we will see, an essentially as good 
fit may be obtained provided A < 0.5. 



Q =8.5 GeV^ 



1.0 - 



0.5 - 



Q =3.5 GeV z 




Q =4.5-5.0 GeV^ 



1.0 
0.5 




Q =2.5 GeV z 



Q =1.5 GeV z 



2.5x10° 10 



Q =0.57-0.58 GeV^ 



0.50 - 



0.25 




2.5x10° 10"' 



g =0.43-0.44 GeV^ 



2.5x10° 10 



Q =0.31-0.34 GeV' 



4x10" 1.6x10° 



4xl0" 4 1.6x10° 



4xl0" 4 1.6x10"- 



Figure 1 
dictions 
HI. Th< 
predictio 
plot is gi 



omparison of pre- 
ith data, Zeus plus 
neutrino data and 
are not shown. The( 
3n of i<2 vs. x 1/l2 . 



We are perfectly aware that we are pushing QCD well below its region of applicability, and that the 
condition of matching at Q 2 — > is at best only of phenomenological value. Nevertheless, the fact that we 
get such reasonable predictions for Ao, po probably indicates that our procedure represents, grosso modo, 
the actual situation, which is also justified by the quality of the fit Eq. (9) provides. If we take all HI and 
Zeus data! 2 ' 3 ! for 0.31 GcV 2 < Q 2 < 8.5 GeV 2 , and we include 10 neutrino xF 3 data 4 we find 

A cS = 0.87 GeV, (e 2 ) B s = 5.28 x 10~ 3 , (e 2 ) B NS = 0.498, (e 2 )C = 0.486, 

for a x 2 /d.o.f. = ^ l 2 4 ■ The value of A e s we have obtained lies somewhere inside the expected bracket, 
A(rif = 2) ~ 0.35 GeV and the value found in the quoted lattice calculation, 0.96 GeV. Clearly, the fit gives 
a compromise, phenomenological quantity. 



4 We choose to include these because the function xFs is pure nonsinglet and its presence thus helps stabilize the 
fits: see ref. 5 for details. 
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The agreement between phenomenology and experiment, shown graphically in Fig.l, is unlikely to 
be trivial; x varies between 6 x 10 and 4 x 10~ 2 , and Fi changes by almost one order of magnitude. To 
see this more clearly, we replace the hard singularity by an evolved soft Pomeron, with a saturated a s . That 
is, we now fit with the expression 



F* = (e 2 q )> 



mog[a s (Ql)/a s (Q 2 )} 



4tt 2 (33- 2n f ) 



+C 




g a s {Q 2 ) 



- d x log 



a s (Q 2 ) 



Q 2 + K s 



+ B NS [a s (Q 2 )]- dNs{1 -^Q 2p s- p 



a s as before. Then we find 

A cS = 0.41 GeV, (e 2 )c = 0.094, (e 2 )C = 0.253, (e 2 q )B NS = 0.41 

and a x 2 /d.o.f. =250/(104-4). We consider these results as convincing proof of the necessity of a hard 
component, which will be henceforth assumed. 

§3. The high Q 2 region. We next discuss the large Q 2 ^ 10 GeV 2 region, under various hypotheses for the 
small Q 2 region, which we then evolve with QCD. We will consider moderately large values of x, x < 0.032 
because we will be interested not only on the leading behaviour as x — > 0, given almost certainly by a hard 
singularity, but also on the subleading corrections. We will, because of the success of the small Q 2 fits just 
discussed, assume for the leading piece a hard Pomeron term with a critical value for A = A . This will be 
also justified a posteriori both by the quality of the fits and by the fact that if, in the best alternative (which 
will turn out to be a hard plus a soft Pomeron), we leave A as a free parameter, its optimum falls very close 
to A ; see below. So we assume that at a certain, fixed Q 2 ~ 1 GeV 2 , one has 

F s (x,Q 2 ) ~ B s x- Xa + F cm .{x,Q 2 ). 

For the correction term, F corr .(a;, Qq), we consider the following possibilities: a soft Pomeron, 



and a P' Regge pole, 5 



F^ IT (x,Qq) ~ constant; 



F^ IT {x,Ql) ~ constant x x 1_ap ' (0) , a P -(0) ~ 0.5. 



We also consider the possibility of adding a triple Pomeron-induced term, 



F TP (x,Qq) ~ constant x x 



-2\„ 



(10a) 
(10b) 

(10c) 



Once assumed the behaviours given in (10), and taking for simplicity that the gluon structure 
function behaves like the quark singlet one, we evolve with QCD for higher Q 2 . Following the methods of 
refs. 5, 6, 8 we find, to NLO, 



F s (x,Q 2 )~B s 1 + 



cs(l + \o)a s 



4tt 



=> gs(l+Ao)a s /47r 



a s (Q 2 



-d+(l+A ), 



Xo +F COII .(x,Q 2 ), (11) 



and the values of the quantities d+, cs, qs may be found in rcf. 5. Then, and depending on the low Qq 
hypothesis we make we find the correction terms, 



h a s (Q 2 ) k a s (Q 2 ) 



Air 



Air 



9£t 



4tt 2 (33- 2n/) 



exp {y/dotr- rfirj , 



(12a) 



For Regge pole theory, and the triple (and multiple) Pomeron terms cf. respectively refs. 13, 14. 
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k, k\ also given in ref. 5. For the P' pole we find a very similar formula: 

F^ T Xx,Ql)~ 1 1 + 2 



4ir 4n 

3 

X \^(i) exp W d °£ t ~ diT ] ■ 



(12b) 



Finally, for the triple Pomeron-induced term, 

F TP (x,Ql) ~ B TP jl + Cs(1 ^ A ° )as }e^ (1+2Ao) "- /4 -[a s (g 2 r d+(1+2Ao) .T- 2A °. (12c) 

As we will see, none of the three possibilities gives a really good fit in the "large" 0.01 < x < 0.032 
region; for the more precise Zeus data the x 2 /d.o.f. is of 1.7. To remedy this we consider the possibility of 
softening the large x region by multiplying Fg(x, Q 2 ) by a factor (1 — x) v : for the discussion and justification 
of this we send to refs. 5, 10. 

The results are summarized in Table I, where we compare the fits obtained with (12a) with the 
fits found in ref. 5 using only the hard singularity or soft Pomeron-dominated expressions, i.e., setting 
respectively cq, Bg equal to zero. 



Table la.- n f = 4; HI data, x < 0.032, Q 2 > 12 GeV 2 



large a; softened | 0.140 GeV 0.318 1.78 x 10~ 4 0.35 



Hard singularity, J A A (e q )Bs (c 2 )Bns xVd-o.f. 

large x softened 

Soft Pomeron only 



110-4 

Ql (e 2 )c {e 2 q )B NS X 2 /d.o.f. 



0.165 GeV 0.70 GeV 0.265 0.246 



110-4 



, j A (fixed) 
1 \ 0.47 1 



?o (e')co (e 2 q )B s (e 2 q )B NS x'/d.o.i. 



Hard + Soft Pomeron* ^ ^ ^ ^ * ^ 



110-4 

Hard + Soft, f A (fixed) Q 2 , (e 2 q )c {e 2 }B s {e 2 )B NS X 2 /d.o.f. 

large x softened* \ 0.47 2.22 GeV 2 0.397 3.36 x 10~ 4 0.353 

* A fixed at 0.230 GeV. 

For the Zeus data one gets similar results, with somewhat less good x 2 /d.o.f. . 

Table 16.- n f = 4; Zeus data, x < 0.025, Q 2 > 12.5 GeV 2 



Hard singularity, j A A (e 2 )Bs (e 2 )Bjvs xVd-o.f. 

only x < 0.01 data \ 0.135 GeV 0.301 1.25 x 10~ 4 0.314 

c u D , / A- Ql (e 2 >c (e 2 q )B NS x7 d -°- f - 

Soft Pomeron only | 0165QeV q gQ Q e y 2 0.282 0.240 ^ 



Hard + Soft Pomeron* 



A (fixed) Ql (e 2 )c (e 2 )B s (e 2 q )B NS x'/d.o.f. 

0.47 2.45 GeV 2 0.252 4.42 x 10~ 4 0.294 



120-4 

Hard + Soft, ( A (fixed) Q% (e 2 )c <e 2 )B s (e 2 q )B NS xVd-o.f 

large x softened* \ 0.47 2.72 GeV 2 0.310 3.417 x 10~ 4 0.370 



A fixed at 0.230 GeV. The optimum value would correspond to A ~ 0.45 GeV. 



In both tables the expression "large x softened" means that we have multiplied the formulas for 
Fs by a factor (1 — x) v , to correct the structure functions for (relatively) large values of x. For the hard 
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singularity case, cf. ref. 5; for the Hard + Soft singularities cases, we have taken v = 10 ~ 11. It is also 
to be remarked that, if we had fitted A to e.g., the Zeus data using a hard plus soft term, we would have 
obtained (not correcting for large x, and fixing A — 0.23 GeV), 

A = 0.429, Ql = 2.40 GeV 2 , (e 2 )c = 0.244, (e 2 q )B s = 3.83 x 10~ 4 , (e 2 q )B NS = 0.30, 

for a x 2 /d.o.f. = 12 ^ 5 . Clearly the values of the parameters, and the x 2 /d.o.f. are so similar to those obtained 
in Table lb fixing A to its low energy value found before, 6 A = 0.47, that we feel justified a posteriori in 
taking such value also at large Q 2 . 

As next possibility, we try a hard Pomeron, plus a P' Regge trajectory, Eq. (12b). If we fit the Zeus 
data with this, adding a softening factor (1 — x) v , v = 11, we find the results of Table II. 



Table II.- n f 


= 4; Zeus data. 






Hard + P' , ( A (fixed) Q 2 , 
large x softened* \ 0.47 1.11 GeV 2 


(e 2 )c P , {el)B s 
0.616 4.25 x 10~ 4 


(e 2 )B NS 
0.343 


X 2 /d.o.f. 

171 

120-4 


* A fixed at 0.230 GeV. 









If we had not corrected for the large x values, i.e., we had not included the factor (1 — x)" , we would have 
obtained a x 2 /d.o.f. of 270. 

We finish this section by describing the results of the fits including a triple Pomeron-induced term. 
That is, we consider that 

F s = F H + F COII .+F TP . (13) 

Then we get, for the Zeus data, 



Table III.- n/ = 4; Zeus data. Triple Pomeron correction added. 



Hard + P J A (fixed) Ql (e 2 q )c P (e 2 q )B s (e 2 )B TP (e 2 )B NS X 2 /d.o.i. 

+TPterm*| 0.47 3.35 GeV 2 0.201 8.67 x 10~ 4 -1.96 x 10" 4 0.315 j±§^ 

* A fixed at 0.230 GeV. 

For HI we would have obtained similar values for the parameters and a x 2 /d.o.f. =105/(110-5). 
The inclusion of the triple Pomeron term improves slightly the fit; but perhaps more interesting is that the 
coefficient of it is, as expected from multi-Pomeron theory,! 14 ! small (when compared with the hard Pomeron) 
and negative. 

Clearly, the best fit is obtained with the hard plus soft Pomerons, just as for the low Q 2 region. 
Not only the x 2 /d.o.f. is quite good, but the values of the parameters are very reasonable. In particular, the 
consistency of the picture is shown by the fact that the value of Q 2 where F^ m is supposed to behave like 
a constant falls in the middle of the low Q 2 region where we showed that precisely this assumptions leads 
to a good fit to data. Because of this, we present in Table IV the parameters of the fits to, simultaneously, 
HI and Zeus data on ep, plus neutrino data. This gives our best set of formulas, providing an excellent fit 
to experiment in a very wide range of Q 2 , x. 



It is of some interest to remark that the optimum A found is actually practically identical to the critical value we 
would have obtained at low Q 2 with zero flavours, A n/=0 = 0.447. The difference with the value for two flavours, 
and the uncertainties inherent to the analysis, however, make us stick to the choice A n ,=2 = 0.470. 
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Table IVo.- n f = 4; Zeus plus HI data; Q 2 > lOGeV 2 , x < 0.01. 



H a rd + P I A ( fixed ) Q « < e ^ Cp < e ^ Bs < e ^ Bjvs ^/d.o.f. 
i 0.47 2.95 GeV 2 0.296 4.28 x 10~ 4 0.349 



144-4 

Hard + P, \ A (fixed) Q 2 , (e 2 ) Cp (e 2 )B s {e 2 q )B TP (e 2 q )B NS X 2 /d.o.f. 

+TPtcrm| 0.47 4.45 GeV 2 0.258 8.33 x 10~ 4 -1.67 x 10" 4 0.359 

A fixed at 0.230 GeV. 

Table IVfc.- n f = 4; Zeus plus HI data; Q 2 > 10 GeV 2 , x < 0.032. x "softened". 



Hard + P 



A (fixed) 
0.47 



2.28 GeV 2 0.311 



(4)B.s 
2.72 x 10" 



(e 2 )B NS xVd-o.f. 
0.315 ^4 



A fixed at 0.230 GeV. 



In the second case (Table IV6) we do not give the fit including a triple Pomeron term as the 
X 2 /d.o.f . docs not vary appreciably if including it provided (e 2 ) \Btp\ 2 x 10~ 4 . We consider the parameters 
given in Table IVa to be the more reliable ones for describing low x structure functions. If we had fitted also 
A with the whole set of data we would have obtained minima for values comprised between 0.42 and 0.49, 
with a variation of the chi-squared of less than two units with respect to the one obtained fixing A = 0.470. 
Finally, if we fit the QCD parameter A, the values which provide minima vary between 0.555 GeV and 0.310 
GeV, and the chi-squared improves by less than five units. Because of this we consider, as stated, that it is 
justified to favour the fits obtained with fixed A = 0.470, A = 0.23 GeV. 
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Figure 2a. Comparison of predictions from Table IV with HI ep data' 2 ' for F2 
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Figure 2b. Comparison of predictions from Table IV with Zeus ep data' 3 ' for F2. 

In Figs. 2 we show the comparison of our fits, with the parameters in Table IV, with data. Note 
that the same values of the parameters are used in Fig. 2a and Fig. 2b. For both Figs. 2 we give the fits 
with the "softened" and straight formulas: the continuous lines indicate large- a; softening, and the dotted 
lines no softening. 

§4. Gluon and longitudinal structure functions. Detailed predictions for the gluon and longitudinal 
structure functions are obtained trivially by adding the soft and hard Pomeron expressions given in ref. 5. 
We will leave the details of this to the reader; but we would like to comment here on a particularly interesting 
prediction of our analysis for the growth of the cross-section <7 7 p_>j/^, p (W) as a function of the cm. energy, 
W . In fact, this cross section may be expressed as a function of the gluon structure function F G , 



a 



M 2 

[ \V) = AF G (x = a-^,Q 2 = M^), 



and A, a are constants approximately known. In the case in which the dominant singularity is a hard one, 
one can obtain, parameter-free the gluon structure function as' 6 ' 

F G (x,Q 2 ) ~ B G [a s ]- d +x- x ° 

x—*0 

and Bq may be calculated in terms of A , B$- So, using above formulas we have, for the logarithmic slope 
of the cross section, 



5 = 



_ logg 7 p^j/y,p(VK) 



2A = 0.94, 



log IT 

The numerical value is obtained with the Ao found in the previous analyses. This is in remarkable agreement 
with the figure reported in a fit! 15 ' including recent HERA data' 16 ' which gives 8 = 0.9. A more detailed 
calculation would produce the whole of F G including corrections sublcading as x — > 0; this unfortunately 
involves a new unknown constant, the soft Pomeron component of the gluon. 



§5. Discussion. The main outcome of our analysis in the present note is that we are able to give a unified, 
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consistent description of small x DIS data, both for large and small values of Q 2 . Besides this, there are a 
number of specific points to which we would like to draw also attention. 

First of all there is the matter of the dependence of our low Q 2 results on the saturation hypothesis 
for a s . It is clear that the good quality of the fits indicates that, with suitable modifications, QCD may 
give a phenomenological description of the data down to very low momenta; but of course this should not be 
construed as a proof of saturation, in particular of the very specific form considered here. One may interpret 
our results, however, as showing that the saturation expression is particularly adapted to represent, in DIS, 
a variety of effects: higher twists, rcnormalons, and likely also genuine saturation. 

A second question is the connection between low and high Q 2 . For the hard piece there is no 
problem, as both expressions are identical up to NLO corrections. For the soft piece, if we start with a 
constant behaviour for Q 2 <~ 2 — 5 GcV 2 , then as Q 2 grows an expression like (2) will start to develop, the 
details of this will depend on what one assumes for the gluon structure function. Because the variation both 
with Q 2 and x of the soft piece is slower than that of the hard part, we think the best procedure is to assume 
constancy of the soft piece up to Q 2 = 8.5 GeV 2 , and the evolved form from there on; since a very good fit 
is obtained at the low momentum region already with the constant behaviour there is little point in adding 
frills, and a new constant (the soft component of the gluon structure function) . 

Finally we devote a few words to the matter of multi-Pomeron exchange. If, at a fixed Q 2 a single 
hard Pomeron gives 7 

F 1P (x,Q 2 ) ~ b 1P (Q 2 >\ 
then an n-Pomeron term will produce the behaviour 

F nP {x,Q 2 ) b nP {Q 2 )s n \ 

s — >oo 

the constants b nP should depend on the momentum at which they are calculated. 

In some approximations (e.g., of eikonal type' 14 ') one has, for Q 2 — — M^ ad , i.e., for on-shell scat- 
tering of hadrons, 

b nP {-ML A ) = (-i) n+1 ^c, 

so we get for the sum 

oo 

F s (x, -M 2 ad ) = £ F nP {x, Q 2 = -M h 2 ad ) 

n=l 

= C- Ccxp \-5s x ] ~ C. 

s^oo 

For values of Q 2 of the order of Qq ~ 2 — 4 GeV 2 , we expect that the b nP (Q 2 ) will not change much, so if 
we write 

b nP {Ql) ~ 6„p(-M^ ad ) + A n . 

we will then get, 

F s (x, Q 2 ) ~ F s (x, -M 2 ad ) + ]T A n ~ C + A lX ~ x + A 2 x~ 2X + ... 

which is the expression we have used in the text. In this respect, the evidence for a behaviour like this may 
be taken as supporting multi Regge theory. 



The following discussion is rather sketchy; details and references may be found in the review of ref. 14. 
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